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Time Averaged Consensus in a Direct Coupled Coherent Quantum 
Observer Network for a Single Qubit Finite Level Quantum System 


Ian R. Petersen 


Abstract — This paper considers the problem of constructing 
a direct coupled quantum observer network for a single qubit 
quantum system. The proposed observer consists of a network 
of quantum harmonic oscillators and it is shown that the 
observer network output converges to a consensus in a time 
averaged sense in which each component of the observer 
estimates a specified output of the quantum plant. An example 
and simulations are included. 

I. Introduction 

There has been significant interest in controlling multi¬ 
agent systems to achieve a consensus; e.g., see [1], [2]. Also, 
the problem of consensus in multi-agent estimation problems 
has been considered; e.g., see [3]. In addition, consensus 
has been considered in quantum multi-agent systems; see 
[4], [5]. The papers [6], [7] considered the problem of 
constructing a direct coupling quantum observer for a given 
quantum system. The problem of constructing an observer for 
a linear quantum system has been considered for example 
in [8]. The theory of linear quantum systems has been of 
considerable interest in recent years; e.g., see [9], [10]. For 
such system models, an important class of control problems 
are coherent quantum feedback control problems; e.g., see 
[9], [11]. In these control problems, both the plant and the 
controller are quantum systems and the controller is designed 
to optimize some performance index. The coherent quantum 
observer problem can be regarded as a special case of the 
coherent quantum feedback control problem in which the 
objective of the observer is to estimate the system variables 
of the quantum plant. The papers [6], [7] considered a direct 
coupling coherent observer problem in which the observer is 
directly coupled to the plant and not coupled via a field as 
in previous papers. This leads the papers [6], [7] to consider 
a notion of time-averaged convergence for the observers. 

We extend the results of [7] to consider a direct coupled 
quantum observer for a single qubit quantum plant, which 
is a network of quantum harmonic oscillators. This quantum 
network is constructed so that each output converges to the 
plant output of interest in a time averaged sense. This is a 
form of time averaged quantum consensus. 

II. Quantum Systems 

Quantum Plant 

We first consider the dynamics of a single qubit spin 
system which will correspond to the quantum plant; see also 
[12], The quantum mechanical behavior of the system is 
described in terms of the system observables which are self- 
adjoint operators on the complex Hilbert space fj p = C 2 . 
The commutator of two scalar operators x and y in Sj p is 


defined as [x : y\ = xy — yx. Also, for a vector of operators 
x in fj p , the commutator of x and a scalar operator y in Sj p 
is the vector of operators [a;, y\ = xy — yx. 

The vector of system variables for the single qubit spin 
system under consideration is 

X p = {X1,X2,X 3 ) T = (<7l, 02, < 73 ), 

where or, a 2 and o 3 are spin operators. Here, x p a vector 
of self-adjoint operators, i.e., x p = x&. In particular x p (0) 
is represented by the Pauli matrices; i.e.. 



The commutation relations for the spin operators are 

[<Jii(Jj] = 2 i^CijkCTk, ( 1 ) 

k 

where denotes the Levi-Civita tensor. The dynamics 
of the system variables x are determined by the system 
Hamiltonian which is a self-adjoint operator on fj p . The 
Hamiltonian is chosen to be linear in x p ; i.e., 

Up = r p x p ( 0) 

where r p G R 3 . The plant model is then given by the 
differential equation 

■Tp(f) — [Xp (f), Up \, 

— A p x p (f), £ p (0) — Xq p , 

z p {t) = CpXp(t) (2) 

where z p denotes the system variable to be estimated by the 
observer and C p £ R lx3 ; e.g., see [12]. Also, A p £ R 3x3 . 
In order to obtain an expression for the matrix A p in terms 
of r p , we define the linear mapping 0 p : C 3 —> C 3x3 as 

/ 0 /3 3 -P 2 \ 

0 p (/3) = I -83 0 A I . (3) 

\ /?2 —/ 3 i 0 / 

Then, it was shown in [12] that 

-i[x p (t),r p x p (t)\ = —2Q p (r p )xp(t) (4) 

and hence A p = —2 © p (r p ). 


In addition, it is shown in [12] that the mapping 0 p (-) has 


the following properties: 

@ P (/3) 7 = -0 P ( (5) 

Q p m 3 = 0 , ( 6 ) 

e P (/3)© P ( 7 ) = 7 P T -P T ii, (7) 

0(©p(/3)7) = 0 p (^)©p(7) - 0 p (7)0 p (/3)- (8) 

Quantum Observer Network The quantum observer net¬ 
work will be a linear quantum system of the form 

x(t) = Ax(t); x(0) = £0 (9) 


where A is a real matrix in R nxn , and x(t) = 

[ X\(t) ... x n (t) ] T is a vector of system observables 

which are self-adjoint operators on an infinite dimensional 
Hilbert space S)\ e.g., see [9]. Here n is assumed to be an 
even number and j is the number of modes in the quantum 
system. 

The initial system variables x(0) = xo are assumed to 
satisfy the commutation relations 

[xj(0),x k (0)] =2iQ jk , j,k = l,...,n, (10) 

where 0 O is a real skew-symmetric matrix with components 
Qjk- The matrix 0 O is assumed to be of the form 

0 O = diag( J, J,..., J) (11) 

where J denotes the real skew-symmetric 2x2 matrix 


The system dynamics (0 are determined by the system 
Hamiltonian which is a self-adjoint operator on the underly¬ 
ing Hilbert space S). For the linear quantum systems under 
consideration, the system Hamiltonian will be a quadratic 
form % = ^x(0) T i?x(0), where R is a real symmetric 
matrix. Then, the corresponding matrix A in <[9]) is given 
by 

A = 2 Q 0 R. (12) 

where 0 o is defined as in (fill) , e.g., see [9]. In this case, the 
system is said to be physically realizable and the commuta¬ 
tion relations hold for all times greater than zero: 

[x 0 (t),x 0 (t) T } = x 0 (t)x 0 (t) T - (x 0 (t)x 0 (t) T ) T 

= 2i0 o for all /. > 0. (13) 

Remark 1: Note that that the Hamiltonian 7~L is preserved 
in time for the system <[9]). Indeed, 'H = ^x T Rx + ^x T Rx = 
— x T RQ 0 Rx+x t RQ 0 Rx = 0 since R is symmetric and 0 O 
is skew-symmetric. 

We now describe the linear quantum system of the form 
© which will correspond to the quantum observer network; 
see also [9], [13], [14]. This system is described by a non- 
commutative differential equation of the form 

x a (t) = A 0 x 0 (t); x o (0) = xoo’, 

z a (t) = C 0 x 0 (t ) (14) 


where the observer output z 0 (t) is the observer network 
estimate vector and A p £ R raoXn °, C 0 £ K : ^ xno . Also, 
x 0 (t) is the vector of self-adjoint non-commutative system 
variables; e.g., see [9]. We assume the observer network or¬ 
der n 0 is an even number with N = ^ being the number of 
elements in the quantum observer network. We also assume 
that the plant variables commute with the observer variables. 
The system dynamics (ITdl) are determined by the observer 
system Hamiltonian which is a self-adjoint operator on the 
underlying Hilbert space for the observer. For the quantum 
observer network under consideration, this Hamiltonian is 
given by a quadratic form: "H 0 = ^x o ( 0 ) T R o x o ( 0 ), where 
R 0 is a real symmetric matrix. Then, the corresponding 
matrix A a in (IT4l) is given by 

A 0 = 2 0 o R o (15) 

where 0 O is defined as in dUV Furthermore, we will assume 
that the quantum observer network has a graph structure and 
is coupled to the quantum plant as illustrated in Figure |T] 



Fig. 1: The graph ( Q,E) for a typical quantum observer 
network. 


The combined plant observer system is described by a 
connected graph (Q, E ) which has N + 1 nodes with node 0 
corresponding to the quantum plant and the remaining nodes, 
labelled 1, 2,..., TV, corresponding to the observer elements. 
This corresponds to an observer Hamiltonian of the form 

7i 0 = ix o (0) T f? o x o (0) 

1 N 

— 7 ^ ^ Xoi ( 0 ) R 0 iX 0 i( 0 ) 

Z z=l 

^ N N 

^ ^ x O j(0) R C ijX 0 j( 0) 

2=1 j =1 

where the vector of observer system variables x 0 is parti¬ 
tioned according to each element of the quantum observer 
network as follows 


x 0 


Xol 

Xo2 


XoN 


We assume that the variables for each element of the quan¬ 
tum observer network commute with the variables of all other 
elements of the quantum observer network; i.e., 

[x 0 2,x 0 j] =0 V i ^ j. 









Here, x ol 


q 0 i 

Poi 


for i = 1,2,... ,N where q 0 i is the 


position operator for the zth observer element and p oi is the 
momentum operator for the ith observer element. 

In addition, we define a coupling Hamiltonian which 
defines the coupling between the quantum plant and the 
quantum observer network: 


N 

H c = x p ( 0 ) T Rc 0 iX O i( 0 ). 
2—1 


Furthermore, we write 


where 


Zol 
Zo 2 

ZoN 


Z 0 i — C 0 iX 0 i ftll i — 1,2, ... , TV . 


Then 


C 0 


Co 2 0 
0 

CoN 


Note that R oi £ R 2x2 , R cij £ R 2x2 , C oi £ R lx2 , and 
each matrix R 0 i is symmetric for i = 1,2,...,TV, j = 
1,2,..., N. In addition, R C Qj £ R 3x2 for j = 1,2 ,... ,N. 
Also, the matrices R C ij for i = 0,1,..., TV, j = 1, 2,..., TV 
are such that R C ij ^ 0 if and only if (i,j) £ E, the set of 
edges for the graph (Q,E). 

The augmented quantum linear system consisting of the 
quantum plant and the quantum observer network is de¬ 
scribed by the total Hamiltonian 

H a = H p + H c + H 0 

1 N 

= r p x p ( 0) + -^x oi (0) T RoiX O i(0) 

i= 1 

1 N N 

+— y ' y ] *02(0) R c ijx 0 j(o) 

2=1 j =1 
N 

+ ^2x p (0) T R c 0 iX oi (0). (16) 

2=1 


Then, it follows that the augmented quantum system is 
described by the equations 


x p {fi) = i[xp (f), Ha \j a^p( 0 ) — XQp‘, 

Xo(l) — ® [Xo (f) > Ha\ j X 0 (C) — Xq 01 

Zp[t) — C p Xp(t ), 

z 0 (t) = C 0 x 0 (t)\ (17) 

e.g., see [12], 

We now formally define the notion of a direct coupled 
linear quantum observer network. 

Definition 1: The matrices R 0 j, R C ij, C 0 j for i = 
0,1,..., TV, j = 1,2,..., TV and the graph ( Q,E) define 


a linear quantum obsen’er network achieving time-averaged 
consensus convergence for the single qubit quantum plant 
(13 if the corresponding augmented linear quantum system 
(fTTb is such that 


lim — 

T-s-oo T 



1 

1 


z p (t) — z 0 (t))dt = 0. 


1 


( 18 ) 


III. Constructing a Direct Coupling Coherent 
Quantum Observer Network 

We now describe the construction of a direct coupled linear 
quantum observer network. In this section, we assume that 
A p = 0 in ©. This corresponds to r p = 0 in the plant 
Hamiltonian. It follows from Q that the vector of plant 
system variables x p (t) will remain fixed if the plant is not 
coupled to the observer network. However, when the plant is 
coupled to the quantum observer network this will no longer 
be the case. We will show that if the quantum observer is 
suitably designed, the plant quantity to be estimated z p (t) 
will remain fixed and the condition <fl8l> will be satisfied. 


We assume that the matrices R C ij, R 0 t for i = 0,1,..., TV, 
j = 1, 2, ..., TV are of the form 

Rcij — U 27 Bj j. Roi — 'Xj / (19) 

where otij £ R 2 , fiij £ R 2 and u>i > 0 for z = 1,2,..., TV, 
j = 1, 2,..., TV. Also, we assume that 

Rcoj = a<yAy where a 0 j = ao =Cp £ R 3 (20) 

for j = 1, 2,..., TV such that (0, j ) & E, the set of edges for 
the graph ( Q,E ). In addition, note that = 0 and /3,j = 0 
for (i,j) E. Furthermore, we assume 

C oi = C p = cxq (21) 


for i = 1,2,..., TV. 

We will show that these assumptions imply that the 
quantity z p (t) = C p x p (t) will be constant for the augmented 
quantum system (fT71 >. Indeed, the total Hamiltonian (H~6l i will 
be given by 

1 N 

H a = -^ 2 uiXoi{ 0 ) T X oi { 0 ) 

Z 2=1 
1 N N 

+ 77 y~~! y~~l Xpi (0) Otij fiijXoj (0) 

2 = 1 j = 1 
N 

+ X p(°) Ta Ojl3ojXoj ( 0 ). 

3 =1 

Now using a similar calculation as in (Q}. we calculate 

Xp(t) = ~i[x p (t),Ha] 

N 

— 2 y ^ 0p(ctQj)tCp(f)/)Q 7 -tr o (f) 

3 =1 

= ~ 2 Q p (a 0 )x p {t) Y $) 3 x o{t)- (22) 

(0j)6 e 










Hence, the quantity z p (t) = C p x p {t ) satisfies the differential 
equation 

z p (t) = -2 C p Q p (a 0 )x p (t) ^2 Poj x <>(t) 

(o ,j)eE 

— @p(o;o)^ : p(f) y ' Pq jX 0 {ty 

(oj)eE 

= 0 (23) 

using ([6| and the fact that 0 ; ,(oj is skew symmetric. That 
is, the quantity z p (t) remains constant and is not affected by 
the coupling to the coherent quantum observer network: 


z p (t) = z p = 2 P (0) Vi > 0. 


(24) 


Also to calculate x 0 (t), we first observe that for any i = 
0,l,...,TV,j = l,2,...,TV. 

\Pij X °j (t) > X oj (t')\ @ij X oj (t)x 0 j (A) X 0 j it) P^jX 0 j ( A ) 

= {Pfj x oj {t)x oj (t) T ) T 

Xoj (A) X 0 j (A ) Pij 
= ( x 0 j(t)x 0 j(t ) ) Pij 

Xoj (t)x 0: , (A) Pij 

= ~ [ x oj (A) , X 0 j (A) ] Pij 

= -2 iJp i:j 

using (fl3T >. Hence, using this result and a similar approach 
to the derivation of (Q2J in [9], we obtain 


Xoj(t) — 


1 [7Aa ) Xoj (A)] 

2 ujjJx 0 j(t) 


1 


N 


T i y ' ( 2 i.JPij) otijXoiit) 


i=i 

N 


+ i ( 2 iJoLji) PjiXoi (A) 


i =1 
„T 


4 ~ iotQjX p (/) ( 2ijPoj) 


N 


2 u>jJx 0 j(t) + J E 

i= 1 
N 

4“ J y ' OtjiPjiXoi^P) 4“ 2 J PojZ p 


(25) 


i=l 


for j = 1,2,... ,7V. 

To construct a suitable quantum observer network, we will 
further assume that 


Xl’ij — , Pij — l^ijOi i (26) 

for A = 1,..., TV, j = 1, 2,..., TV where ( i , j) £ E. Here, 
ay £ M 2 and 

/^ i] = Liji T> 0. (27) 

Also, we will assume that 

Ay = — AHyou (28) 

for j = 1, 2,..., TV where (0, j) £ £7. 


In order to construct suitable values for the quantities /Mj 
and uji so that (fl8l) is satisfied, we will require that 

2cz j J a ]_ ^a i a i a i 

(i,j)£E,i>0 

— ^2 fiij Ja\a{ai 4- 2JpojOt{ai = 0 (29) 

(i,j)£E,i>0 

for j = 1, 2,..., TV. This condition is equivalent to 

Uj = ^2 Myl|ai|| 2 + Moil|cti|| 2 (30) 

(i,j)EE,i> 0 


for (0, j) £ E and 

LOj = ^2 My II ay 

(i,j)£E,i> 0 


(31) 


for (0 ,j)?E. 

Then, we define 

Xoj {t) = Xoj(t) — 1 | 2 CXlZp 

for j = 1,2,..., AT. It follows from (|29t and (l25t that 

N 

— 2ujjJx oj (t) H - J ^ ^ fiij Q'i j Xpj (t ) 

i=l 

N 

T J ^ ^ OLji^j^Xoiij^) 

i— 1 

— e ZUjJx 0 j{t) 2 ^ ^ HijJoixOLx X 0 i (f) 

( 2 5 i 7 ')G-E,i >0 


for j = 1,2, 

We now write this equation as 


’ Xoi {t) 


Xol (t) 

Xo2(t) 

= A 0 

Xo2 (t) 

_ X oN {t) 


_ X 0 N{t) 


where A 0 is an iV x TV block matrix with blocks 


(32) 


O'oij 


2 ujiJ for i = j, 

—2^ijJa\aJ for i ^ j and (i, j) £ E, 

0 otherwise 


for i = 1, 2,..., TV, j = 1, 2,..., TV. Also, A a is as given 
in (Q2J where R a is a symmetric TV x TV block matrix with 
blocks 


uy/ 


Tnii = < — 


Mycyaq 


for i = j, 

for -t ^ j and (i, j) £ E, 

otherwise 


for i = 1,2,. ,.,TV, j = 1,2, ...,TV. 

To show that the above candidate quantum observer net¬ 
work leads to the satisfaction of the condition (fl8l) . we note 


that 


x 0 


Xol 

Xo2 


XoN 









where 


satisfies (|32| >. Hence, if we can show that 


1 ( T 

Inn - / x 0 (t)dt = 0 

T—>-°o J Jq 


(33) 


11*01 || 
ll*o 2 || 


then it will follow from 


|| x oN 


c n 



a l 

ai 


0.1 



s 


01 

1 

aj 0 


01 

IMI 2 

0 




_ T 




a i . 


L o 1 j 



(34) 


and R 0 is a symmetric TV x TV matrix with elements defined 
by 

( tOi for i = j, 

-faj for i j and (i.j) £ E, 

0 otherwise 

for i = 1,2, j = 1,2, 

Now the vector x 0 will be non-zero if and only if the 
vector x 0 is non-zero. Hence, the matrix R 0 will be positive- 
definite if we can show that the matrix R 0 is positive-definite. 
In order to establish this fact, we first note that (l30l > and (OH 
imply that 

U}j = '^2 dij + do j 


that m is satisfied. 

We now show that the symmetric matrix R a is positive- 
definite. 

Lemma 1: The matrix R a is positive definite. 

Proof: In order to establish this lemma, let 


for (0, j) £ E and 

= E dij\\oi\\ 2 

(i,j)eE,i> o 

for (0, j) ^ E. Hence, we can write 


x 0 


Xol 

Xo2 


XoN 


be a non-zero real vector. Then 


X~ R n Xa 


> 


N 

J 2 ^\\Xoif 

E dij x oi a l x oj a l 
(i,j)£E,i>0,j>0 
N 

J2^W X oi\\ 2 


E wjii Xoi\\\\x 0 j\\ II Oil I 

( i,j)£E,i>0,j>0 


N 

E 

2=1 


Ldi\\X 0 


^ ^ ll^oz || \\ x oj | 

(i,j)£E,i>0,j>0 


(35) 


using the Cauchy-Schwarz inequality. Here, 

dij = MMI 2 

for 0 = 1,2,..., N, j = 1,2,..., TV. Thus. d35l > implies 


x 0 RqX 0 ■ 1 x Q R 0 x o 


Ro — Rol "F Ro2 


where R 0 1 is a symmetric TV x TV matrix with elements 
defined by 


^ olij — 


Ei(k,j)eE,k >0 dkj 
~dij 

0 


for i = j, 

for i f j and (i.j) £ E, 
otherwise 


for i = 1, 2,..., TV, j = 1, 2,..., TV. Also, R a 2 is a diagonal 
TV x TV matrix with elements defined by 

_ ( doj for i = j and (0,j) £ E, 

° 2 ^ \ 0 otherwise 


It follows that the matrix R a 2 is positive semidefinite. 

Now the matrix R 0 1 is the Laplacian matrix for the 
weighted graph {G,E) obtained by removing node 0 from 
the graph (Q, E) along with the associated edges. Then each 
edge (i. j) £ E is given a weight fij J e.g., see Figure 
U which shows the weighted graph (Q, E) which would 
correspond to the graph (Q, E) shown in Figure Q] 

It follows that the matrix R a \ is positive-semidefinite with 
null space of the following form: 

N(Roi) = span{/i, / 2 ,..., fm} 


where m is the number of connected components of the 
graph (G,E). Also, each of the vectors fi, f2, ■ ■ ■, fm are 
vectors whose elements are either zeros or ones. For the 
vector fk, the elements of this vector which are ones corre¬ 
spond to the nodes in the graph (Q 1 E) in the fcth connected 
component. 
















1 


5 


10 


for all x o (0) and t > 0. Hence, 



Fig. 2: The weighted graph (G,E) corresponding to the 
graph (Q, E) in Figure [0 


The fact that R 0 \ > 0 and R a 2 > 0 implies that R 0 > 0. 
In order to show that R 0 > 0, suppose that x is a non-zero 
vector in Af( R a ) ■ It follows that 

x T R a x = x T R 0 \X + x T R a ix = 0. 

Since R 0 \ > 0 and R 0 2 > 0, x must be contained in the null 
space of R ol and the null space of R 0 2 . Therefore x must 
be of the form 

m 

fc=1 

where not all 77 = 0. However, since the graph (G, E) is 
connected, it follows that there must be at least one branch 
(0. j) £ E to a node in each of the connected components 
in the graph (G,E). Then 

X T Ro 2X = k{j) = 0 

(0 ,j)GE 

where k (J ) corresponds to the node of the connected com¬ 
ponent in ( G,E) which the branch (0 ,j) connects to. Since 
each p, a j > 0 , it follows that 

7 k(j) = 0 

for all (0 ,j) € E. Furthermore, since each connected 
component in ( G,E ) has at least one branch (0,j) G E 
connected to it, it follows that 71 = 72 ■ ■ ■ = 7 m = 0 . 
However, this contradicts the assumption that not all 77 = 0. 
Thus, we can conclude that the matrix R a is positive definite 
and hence, the matrix R a is positive definite. This completes 
the proof of the lemma. ■ 

We now verify that the condition ( fl 8 l > is satisfied for the 
quantum observer network under consideration. We recall 
from Remark Q] that the quantity ^x 0 (t) T R 0 x 0 {t) remains 
constant in time for the linear system: 

x 0 — -A 0 x 0 — 2QR 0 x 0 . 

That is 

^x 0 (t) T R 0 x 0 (t) = ii o (0) T f? o i o (0) Vf > 0. (36) 

However, x 0 (t) = e 2 eRot x o (0) and R a > 0. Therefore, it 
follows from (l36t that 

V / A mm (f? o )||e 2efl ° t x o (0)|| < VA mM (i? o )||i o (0)|| 


ll e 


2QR 0 t\ 


< 


A max (Ro) 

^min (Ro) 


(37) 


for all t > 0. 

Now since 0 and R 0 are non-singular, 

rT 


/o 


e 2BR 0 t dt = l e 2eR 0 T R -l Q -l _ 

2 ° 2 ° 


and therefore, it follows from ( 177 b that 

rT 




e 29Rat dt\\ 


h\l e 2@R 0 T R -lQ-l _ l R - 10 - 1 | 
T 2 ° 2 ° 1 


< 


< 


2QR a T | 


2T 


I R-'Q- 


2 T 


IK -1 © 


’ ^max (Ro ) „, D —— 11 


1 


2 T V A min(Ro) 


I \RoO~ 


+^\\Ro 1 o ~ 1 1 


-»• 0 


as T — > oo. Hence, 

lim 4II / x 0 (t)dt\\ 

t-> oo 1 J 0 

= 1 ™ ^11 [ T e^xMdtW 

T—too T J o 

< lim i|| [ T e 2eRot dt\\\\x o m 

oo T J o 

= 0 . 


This implies 


1 ( T 

hm - / Xo (t)dt = 0 
r-»oo T J o 


and hence, it follows from (l32t and ( 1341 that 


lim — 
T- s-oo T 



Also, (l24t implies 


lim - / z p (t)dt 

T->oo 1 J Q 


1 

1 

1 


1 

1 



Therefore, condition (ITRl i is satisfied. Thus, we have estab¬ 
lished the following theorem. 

Theorem 1: Consider a single qubit quantum plant of the 
form (0 where r p = 0 and hence A p = 0. Then the 
matrices R ci , R C ij, C 0 i, R 0 i for i = 1,2, ...,7V, j = 
l,2,...,iV and the connected graph (G,E) will define a 
















direct coupled quantum observer network achieving time- 
averaged consensus convergence for this quantum plant if 

the conditions ED, ED- GD. (ED. ED, ED- ED- ED 

are satisfied. 

Remark 2: The quantum observer network constructed 
above is determined by the choice of the positive parameters 
fj,ij for i = 0,1,..., AT, j = 1,2, ...,7V. A number of 
possible choices for these parameters could be considered. 
One choice is to choose all of these parameters to be the 
same as jlij = wo for i = 0,1,..., TV, j = 1, 2,..., TV 
where wo > 0 is a frequency parameter. Another possible 
approach is to choose the parameters for i = 0,1,..., TV, 
j = 1,2,..., TV randomly with a uniform distribution on a 
suitable frequency interval. 


can be described by the equations 


x a (t) = A a x a (t), where x a (t) 


Zp(t) 

X 0 l{t) 

Xo2(t) 


and 


*^o5 (t) 


oo o 
0 0 10 

2 -10 0 

0 0 0 

2 2 0 

0 0 0 

2 2 0 

0 0 0 

2 2 0 

0 0 0 

2 2 0 


0 0 0 

0 0 0 

2 0 2 

0 10 0 

-10 0 2 

0 0 0 

2 0 -10 

0 0 0 

2 0 2 

0 0 0 

2 0 2 


0 

0 

0 

0 

0 

10 

0 

0 

0 

0 

0 


0 0 

0 0 

2 0 

0 0 

2 0 

0 0 

2 0 

0 10 

-10 0 

0 0 

2 0 


0 0 

0 0 

2 0 

0 0 

2 0 

0 0 

2 0 

0 0 

2 0 

0 10 

-10 0 


IV. Illustrative Example 


Then, we can write 

X a (t) = $(t)x a ( 0) 


We now present some numerical simulations to illustrate 
the direct coupled quantum observer network described in the 
previous section. We choose the quantum plant to have A p = 
0 and C p = [ 1 0 0]. That is, the variable to be estimated by 
the quantum observer is the spin operator cri of the quantum 
plant. For the quantum observer network, we choose TV = 
5 so that the quantum observer network has five elements. 
Also, we suppose that the graph (Q, E) defining the plant 
observer network is the complete graph corresponding to the 
five observer nodes and the plant node; i.e., every node is 
connected to every other node in this graph. This graph is 
illustrated in Figure [3] In addition, we choose a± = [1 0] T 
and as discussed in RemarkEl we choose the parameters fi, :i 
so that jlij = u>o for i = 0,1,..., TV, j = 1,2,..., TV where 
wo = 1. Then the dynamics of the corresponding quantum 
observer network are defined by equations ED and ED . 



Fig. 3: The plant observer network considered in the exam¬ 
ple. 


where 

$(f) = e Aat . 

Thus, the plant variable to be estimated z p (t) is given by 
z p {t) = e\C a <S>{t)x a (A) 

2N+2 

= ^2 eiC a ®i(t)x ai (0) 

i= 1 

where 



e\ is the first unit vector in the standard basis for R w+1 , 
<!>i(f) is the 7 th column of the matrix <f>(f) and x Q i(0) is 
the 7th component of the vector x„.(0). We plot each of 
the quantities e 1 C a $ 1 {t),e 1 C a $ 2 {t), ■ ■ ■, e 1 C a $ 2 N+i(t) in 
Figure [4ja). 

From this figure, we can see that eiC a $i(f) = 1 and 
eiC a &2{t) = 0, eiC' a $2(f) =0, ..., eiC a <&2N+2(t) = 0, 
and z p (t) will remain constant at z p (0) for all t > 0. 

We now consider the output variables of the quantum 
observer network z Q i(t) for i = 1,2,..., TV which are given 
by 

2 AT +1 

Zoi (f) — ^ ^ doA^j {t)Xaj (0) 
i =i 

where a+i is the (i + l)th unit vector in the stan¬ 
dard basis for M Ar+1 . We plot each of the quan¬ 
tities e i+1 C a $ 1 (t),ei +1 Ca$ 2 (t),...,e i+1 Ca$ 2 N+ 2 (t) in 
Figures Hb) - fflT). 

To illustrate the time average convergence prop¬ 
erty of the quantum observer ED, we now plot the 
quantities A e i+1 C a ®i(t)dt, A f Q e i+1 C a $ 2 (t)dt, ..., 
A f2 e l+1 C a $ 2 N+ 2 {t)dt for i = 1,2,..., TV in Figures 
|Da)EI e )- These quantities determine the averaged value of 
the ith observer output 

. -T 2N+1 

z a o r(T) = - / ei+iC a ^(t)x aj (0)dt 

1 Jo 3=1 


For this example, the augmented plant-observer system 
















Fig. 4: Coefficients defining (a) z p (t), (b) z 0 \{t), (c) z 0 2 {t), 
(d) z o3 (t), (e) Zoi{t), and (f) z o5 (t). 


for i = 1,2 ,,N. From these figures, we can see that for 
each i = 1,2 ,N, the time average of z 0 i(t) converges 
to z p (0) as t —>• oo. That is, the quantum observer network 
reaches a time averaged consensus corresponding to the 
output of the quantum plant which is to be estimated. 
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